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Abstract 

A new set of hadronic form factors, which has been implemented in TAUOLA, is described. 
Keywords: Hadronic decays of the Tau lepton, Monte Carlo methods. Resonances, Chiral symmetry. 



1. Semileptonic tau decays in TAUOLA 

Tau decays including hadrons are a privileged sce- 
nario HI] to study non-perturbative QCD in a rather 
clean environment provided by the electroweak half 
of the process. In Ref. [2] it was concluded that the 
most essential missing step to perform, concerning the 
study of hadronic tau decays, was the appropriate choice 
of the hadronic currents. The original version of the 
Monte Carlo (MC) generator TAUOLA H used the so- 
called Kiihn-Santamaria (KS) model f^, and its exten- 
sions, to construct them. Within this model, they are 
built to fulfil the leading order (LO) result in the low- 
energy effective field theory of QCD, Chiral Perturba- 
tion Theory (xPT) [5], but they violate the next-to- 
leading order (NLO) one 0]. This approach was suf- 
ficient and succesful twenty years ago, but already the 
CLEO and Aleph Collaborations realized, later on, de- 
partures of the predictions from data, a feature which 
could be expected taking into account several incon- 
sistencies in later parametrizations of the three meson 
modes including Kaons This resulted in private 
versions of the code, with fine-tuned initializations - 
which sometimes violated basic principles of QCD|@]- 
that were documented in Ref. |@] . Nowadays, with the 
massively increased data samples from the B-factories 
BaBar and Belle -the most of which have not been an- 
alyzed yet-, it is pressing to upgrade the hadronic cur- 
rents in TAUOLA in order to obtain as much QCD in- 
formation as possible from experiment; moreover with 
the perspective of the super-flavour factories producing 



huge amounts of high-quality data in the near future. 

In Ref. 1 10] a new set of form factors for hadronic tau 
decays based in analytical results obtained from Reso- 
nance Chiral Theory (RvT nl ill is documented. In these 
Proceedings, a description of the implementation of the 
new modules of the MC program for its user is given 
in Ref. lll2ll . Here we focus on the hadronic currents 
themselves. 

Lorentz invariance determines the most general de- 
composition of the hadronic current. This is 

r = N[{iH - PiYF^is) + (pi + pifF'is)] , (1) 

in the two-meson channels, where F^(s) and F^{s) are 
the vector and scalar form factors, given in terms of s = 
(Pi + Pif, and 

r = N{n[(p2 - PiYFi + (p3 - PiYFi (2] 



+(pi - piYFi] + q^F4 ■ 



7^vpa-P\P2P3F5 



in the three-meson decayfl where T^v - gfiv - q^qvlq^ 
stands for the transverse projector, and - {pi + pi + 
PiY is the momentum of the hadronic system. Only 
two among the Fi, F2 and Fj, (axial-vector) form fac- 
tors are independent. F4 is the generally suppressed 
pseudoscalar form factor and F5 is the vector form fac- 
tor The kinematical invariants are q^, and two of the 



The (by-far) most important four-meson tau decay, into four pi- 
ons, is currently parametrized following Ref. fT3h . 
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three = (/?,■ + pj)^ - s^, i + j k - 1,2,1), since 
S3 - — s\ — $2 + m\ + + m^. Chiral symmetry 
relates the N factors for the diverse decay channels in 
Eqs.([T]) and (|2|i, respectively. Unfortunately, there is no 
analytic way of deriving the expression of these form 
factors from the QCD Lagrangian. This does not mean 
that the underlying theory is useless to find them, as we 
discuss next. 



2. Theoretical setting 

QCD has a well-defined expansion parameter at low 
energies within the light quark sector, h^pj, that allows 
to build ;t'Pr. The associated approximate chiral sym- 
metry is useful to understand the data at low values of 
the invariant mass of the hadronic system, but it is in- 
sufficient to explain them throughout the whole phase 
space [14]. This happens because A.^pt is no longer 
small for E > nix- 1 /Nc 1, 15 1 seems a good candidate to 
build the expansion upon, given its success in explaining 
many features of meson phenomenology Further- 
more, it proves efficient in dealing with effective field 
theories of QCD in the low- and intermediate-energy 
regions, ^PT and RxT in our context fvf]. Working 
at lowest order in 1/Nc amounts to consider an infi- 
nite number of stable resonances which experience lo- 
cal effective interactions at tree level among them. This 
supports RxT, which includes the xPT Lagrangian at 
LO, and reproduces the one at NLO upon integration of 
the resonances [Tl'l, which are active degrees of free- 
dom in the theory. R^T provides a rigorous way lIlSIl 
of computing a NLO effect in this expansion, the reso- 
nance widths, of paramount importance to understand 
hadronic tau decay data. A model independent real- 
ization of the infinite tower of resonances remains un- 
known. We model this setting by cutting the spec- 
trum in a way that resembles nature. One hopes that 
since the lowest-lying resonances dominate the light- 
flavour Physics, this procedure would allow to capture 
the essentials of the involved dynamics. The RxT, de- 
rived from symmetries, is still lacking the short-distance 
QCD behavior. When it is imposed to the Green func- 
tions 1 19] and associated form factors, a number of re- 
lations between the Lagrangian couplings arise, making 
the theory more predictive. 



s « Mp. Computations at NNLO are available both in 
the S U(2) 120] and in the S U(3) fll] symmetry cases. 



In the region Mp 



< s < 



1 GeV, the chiral expansion 



breaks down and the dominant p(770) exchange has to 
be accounted for Several approaches have been de- 
veloped. Among them, matching xPT results to vec- 
tor meson dominance using an Omnes solution ll22ll for 
the dispersion relation [23], employing an Omnes so- 
lution for the dispersion relation |24] or utilizing the 
unitarization approach 0251]. For larger energies, in the 
1-2 GeV region, the excited resonances play an impor- 
tant role and shall be incorporated to the description. 
Ref. i6] includes the p(1450) through a Schwinger- 
Dyson-like resummation and Refs.[27] include a tower 
of resonances inspired from dual QCD. Since Ref. 1123 1 
will be our starting point, let us recall their main features 
in the following. 

In this case - V2 in Eq.([T]), and the scalar form 
factor is zero in the S U (2) symmetry limitH. The vector 
form factor, at NLO in xPT, is 



F\s) = 1 + 



-s- 



Fl 96nFl 



A„{s) + ]^Ak{s) 



with 



Apis) = Log 



+ cTpLoa 



, (3) 



{ap+V 
Wp-I, 



a-p = A 1 - 



4m2 



(4) 



The computation in RxT, within the antisymmetric ten- 
sor formalism, reads 



F^s) = 1 + ^^^^ ' 



(5) 



When F (s) > is required, the condition FyGy — 

FI is found, which yields the vector meson dominance 
prediction (/i ~ Mp) 



Ml 
Mj-s 



IMj 



(6) 



The matching of Eqs.O]) and ^ is straightforward 



Ml 



Mj - s 96nFl 



An{s) + -Ak{s) 



(7) 



There are different approaches to deal with the di- 
verse energy regimes which are probed through this de- 
cay: xPT should be a valid description of the data for 



-Even when first order isospin violating corrections are included 
12811 it does not contribute. 
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When unitarity and analiticity properties are required, 
the Omnes solution emerges 



Ml 



Mj- s 



exp<- 



967iFl 



(8) 

Not surprisingly, the p{110) off-shell width is related to 
the imaginary part of the same loop function 



-MpS 



AAs) + ]^Ak{s) 



(9) 



A possible solution to avoid double counting of the 
imaginary parts, which was adopted in Ref. [23J, is 



Mi 



Ml-s- iMpVpis) 



r'Rc[A,(s)+|Ax(s)]} 



(10) 

Eg. ([Toll reproduces xPT at NLO, vanishes at s ^ 
cx), has SU(2) symmetry built-in and complies with 
analiticity and unitarity constraints up to first order in 
the expansion of the exponential. This description was 
successfully confronted to data using only one parame- 
ter, Mp. Present data have become much more precise 
and the Belle results 129*1 point to an interference pattern 
between excited resonances in this decay. All this moti- 
vates us to include Isoll . analogously, the contribution of 
the excited resonances [p' = p(1450) and p" - p(1700) 
in this case] while keeping these nice properties Bill 



M2 



s - iMpTpis) 



' Ml-s- iMp.T pis) 



5se"i''- 
Ml, - s-iM„nTp„{s) 



The parameters y and 6 are related to RyT couplings for 
the excited resonances [as yKn in Eq.(fT3ll1. the Tp'{s) 
and rp"(i) widths are modeled as decays to two pi- 
ons, and the phases (p\ and <p2 should vanish, at least, 
as 1/A^c- Eq. lfTTl i corresponds to what is included in 
TAUOLA right now 0. SU(2) breaking has only been 
coded partially, through the kinematical and loop func- 
tions. Electromagnetic corrections ||28, 32] have been 
considered llsoll but not incorporated to the MC yet. 



^^The values of the parameters for this mode and the others can be 
found in the quoted references, where comparisons to data are also 
available. 



An alternative approach to using the Omnes solution 
consists in employing an n-subtracted dispersion rela- 
tion where the relevant phaseshift, 6\{s), is obtained as 
V)mF^ {s) I'^KeF^ (s). In this procedure llssll . unitarity 
and analiticity are satisfied to all orders with ifsoll 



F\s) = 



Ml + i(ye''*' +5e'^^-) 



Ml[l+ ^p%e [a,{s) + {Ads))] - s - iMpFpis) 
yse'"!"' 

Mj, [l + ^p.%eAAs)\ - s - iMp.Tp,{s) 
6se'<l>'- 

Ml„ [\ + ^p.,%eAAs)\ - s - iMp,.Yp.,{s) 



(12) 



in which 



and ^p 



KM\crl(M\ 



(analogously 



for ^p"). The result for the resummation in Ref. OlSll has 



been employed in the denominator of the p contribution. 



4. Other two meson t decay channels 

The — > K^K'^Vj decays are again described only 
in terms of the vector form factor to an excellent de- 
gree of approximation. The current parametrization 
in TAUOLA follows the Guerrero-Pich formula, see 
Eq.lfTOb 123 . 34]. There is also an option to use Eq. lfTTT i. 
Further developments in F^^{s) will be immediately 
translated to F^^{s). The vector form factor in the 



Ref. Q 



(Kn) Vt decays is currently coded following 



Ml. + SyKn 



SyKn 



Ml,-s- iMK-FK-is) Ml„ - s - iMK-'TK-'is) j 



exp 



UStt^F^ 



%e[AKjs)+AK,(s)] 



(13) 



^^TTie function A/)q(5) is ^ 



Apq(s) 



192n^sMpQ(s)-LpQ{s)] 



(14) 



in the notation of Gasser and Leutwyler. An option will 
be given to switch between this form factor and the one 
in Ref. ifssll Fl. The scalar form factor is fundamental in 
this decay channel to achieve a precise description of the 
decay data fyi^ at low values of s. Moreover, it is es- 
sential to understand CP violation in this channel, which 
has been reported recently jssll . TAUOLA is ready ^ 
to handle such and distinguishing terms. The im- 
plementation of F^ (s) [,39,1 is documented in Ref. ll40ll . 



Other interesting approaches are those of Refs.l36 
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5. Three meson r decay channels 



The 



(nnn) Vj and t 



(KKtt) Vr decays 



have been coded following Refs. 114 111 . One- and two- 
resonance exchange diagrams were considered within 
RxT and the appropriate short-distance behaviour was 
required, yielding sets of compatible relations among 
the Lagrangian couplings in both decays (including 
FyGy = F^, as in the two meson tau decays). The 
p' resonance was introduced phenomenologically to im- 
prove the description of the data in the {mmYvr 
decays. The progress with respect to the earlier descrip- 
tion given by the KS model can be appreciated in Fig. 
1 of Ref . f42] . The inclusion of final state interactions 
(FSI) in these decays is under study |40]. It should im- 
prove the agreement with data in the dF/dsij distribu- 
tions, specially at low values of sij. 

6. Conclusions 

A set of form factors based on RxT calculations, 
corresponding to 88% of the hadronic width of the r 
lepton, has been implemented in TAUOLA. They are 
ready for precise confrontation with data gathered at 
Belle and BaBar (and future Belle II & Frascati superB 
facilities). In order to obtain the maximum possible 
information from experiments, the theory input to 
the MC has to be as accurate as possible with known 
properties respected (xPT results at low energies, 
smooth behaviour of the form factors at short distances, 
unitarity, analiticity, ...). Still, there are improvements 
to be done in all modes: appropriate inclusion of SU(2) 
breaking in the tt^tt'-'' channel, stabilization of F^^(s), 
inclusion of excited resonances in the KKn modes, 
and addition of FSI (mainly the cr effect) in the Stt mode. 

I congratulate the organizers of PHI PSI 11 and 
of the WG on MC Generators for low-energy Physics 
for their excellent job. I acknowledge funding of the 
FPA2007-60323 and CPAN (CSD2007-00042) grants. 
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